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ABSTRACT 

Thermally pulsing asymptotic giant branch models of globular cluster stars are cal- 
culated using a synthetic model with the goal of reproducing the chemical compo- 
sition, core masses and other observational parameters of the four known globular 
cluster planetary nebulae as well as roughly matching the overall cluster properties. 
The evolution of stars with an enhanced helium abundance (Y) and blue stragglers are 
modeled. New pre-thermally pulsing asymptotic giant branch mass-losses for red giant 
branch and early as ymptotic g iant branch stars are calculated from the Padova stellar 
evolution models (jBertelli et all [20081 120091) . The new mass- losses are calculated to 
get the relative differences in mass- losses due to enhanced helium abundances. 

The global properties of the globular cluster planetary nebula are reproduced with 
these models. The metallicity, mass of the central star, overall metallicities, helium 
abundance and the nebular mass are matched to the observational values. Globular 
cluster planetary nebulae JaFu 1 and JaFu 2 are reproduced by assuming progenitor 
stars with masses near the typical main sequence turn-offs of globular clusters and 
with enhanced helium abundances very similar to the enhancements inferred from fit- 
ting isochrones to globular cluster colour-magnitude diagrams. The globular cluster 
PN GJJC-1 can be roughly fit by a progenitor star with very extreme helium enhance- 
ment {Y PS 0.40) near the turn-off producing a central star with the same mass as 
inferred by observations and a very low nebular mass. The abundances and core mass 
of planetary nebula Ps 1 and its central star (K648) are reproduced by a blue straggler 
model. However, it turned out to be impossible to reproduce its nebular mass and it 
is concluded some kind of binary scenario may be needed to explain K648. 

Key words: stars:AGB and post- AGE - blue stragglers - mass-loss - white dwarfs 
- globular clusters: general - planetary nebulae: individual 



1 INTRODUCTION 

The globular cluster (GC) system of the Galaxy contains 
four known planetary nebulae (PNe) all of which have un- 
usual abundances. The nebulae Ps 1 and its central star 
K648 in M15 ha v e a very high C/O ratio of nearly 10 
llAlves et al.ll2000l: iBianchi et all 1 19951. l200ll : [Howard et all 
[ 19971 : [Ranch et al.ll2002f ). The PNe JaFu 1 in Pal 6 and JaFu 

2 in NGC 6441 have high ratios of He/H (|jacoby et a l.| |l997| ) 
(0.115 and 0.141 for JaFu 1 and JaFu 2, respectively). The 
fourth GC PN GJJC-1 in M 22 is one of three kn own PN with 
no hydrogen in the nebula ([Jacobv et al.[[l997f ). All four are 
very different from the typical disc PN. They are even more 
difficult to explain since the typical age of a globular clus- 
ter is ~ 12.5 Gyr leading to a turn-off mass of ~ 0.85 Mq. 



E-mail: buelljf@alfredstate.edu 



At this mass only minor changes from the zero-age main 
sequence (ZAMS) abundances to PN abundances are ex- 
pected. If the progenitors had a typical halo or thick disc 
compositon the abundance ratios would be similar to those 
found in disk PNe (He/H~ 0.100, C/0< 1 and N/0~ 0.4). 
The origin of these PNe requires an explanation. 

[Jacobv et ah I l|l997h surveyed 130 GCs in a search for 
PN in GC. They concluded that the number known is prob- 
ably complete. They found this is less than the expected 
number (« 16) given the total luminosity of the Galactic 
globular cluster system. This number is based on stellar 
populations which are younger than the stars in globular 
clusters. GCs have very low mass turn-offs (~ 0.85 Mq). At 
this mass a typical PN central star (CSPN) should have a 
mass of ~ 0.52 Mq. The time required from when a star of 
this mass leaves the asymptotic giant branch (AGB) for this 
star to get hot enough to ionize the ejected mass would be 
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too long. The ejected mass would dissipate before a visible 
PN cou ld be observed . This is known as a "lazy" planetary 
nebula (jRenzinill 19791 ). 

The masses of the central stars of the globular clus- 
ter PN are larger than the typical measured masses of 
the youngest white dwarfs (WDs) at the top of the 
cooling sequence. The cooling sequences of white dwarfs 



the expectations of the PN phase from all the generations 
of GC stars. I also model the expected PN phase of the blue 
straggler stars. In section 2 I describe the TP-AGB model 
used. In section 3 I describe the results of these models. In 
section 4 I discuss some implications. In section 5 I summa- 
rize the results. 
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ity of main sequence stars evolve into white dwarfs of mass 
0.5-0.55 Mq. and the average of a young GC white dwarf is 
probably between 0.50 and 0.53 Mq. AU of the G C PNe have 
central star (CSPN) mass e s abo ve this average. lAlves et al.l 
l|2000l) and iBianchi et all (|200ll ') find t hat the mass of the 
central star of Ps 1 is 0.58-0.60 Mq. jjacobv et al.lll997l ') 
determined the CSPN of JaFu 1 and JaFu 2 have masses 
of 0. 55 MfT). The estim ated mass of GJJC-1 is around 0.56 
Mq jPefia et al.lll992l ). 

Both the CSPN masses higher than the white dwarfs 
and the unusual abundances seem to require unusual stars 
as progenitors. GCs contain a variety of unusual star types 
in addition to the standard types of stars which should be 
considered as potential progenitors of GC PNe. GCs contain 
blue stragglers which are thought to be the mergers of two 
main sequence stars and hence act like main sequence stars 
with masse s higher then the turn-off mass. iBianchi et al.l 
l|200ll ) and lAlves et al.l (|2000t ) suggested the projenitor of 
Ps 1 is a blue straggler because of its high core mass and 
evidence of at least one third dredge up (TDU) event as 
evidenced by the very high C/O ratio in the nebula. The- 
oretical models of thermally pulsing- AGB (TP-AGB) stars 
suggest a core mass of ~ 0.58 Mq is required to get a TDU 
event. 

Another type of star found in globular clusters are sec- 
ond generation stars which incorporate material from the 
more massive stars of the first generation (primordial com- 
ponent or P). The first generation stars have abundances 
which reflect the abundances of the interstellar medium 
material from which they formed and hence have a nor- 
mal amount of helium {Y « 0.25). The second genera- 
tion often incorporates material with a higher fraction of 
helium than normal {Y ~ 0.30). Some clusters show evi- 
dence for additional populat i ons w ith even higher Y val- 
ues (E.g. ICaloi fc D'Antonal (|2007h '). These multiple pop- 
ulations show up observationally in a number of ways. It 
shows up as an Na-O anticorrelation in both red giants 
and in main s equen ce stars. This was flrst observed by 
iGratton et al.l (|200lh when they noted the Na-O anticorre- 
lation shows up in main-sequence stars in addition to red gi- 
ants in several clusters. Some clusters have distinct multiple 
main sequences (e .g. ujCen (Bcdin ct al. 2004), NGC 2808 
jPiotto et al.ll2007l ). etc.) which can be flt by having a second 
main sequence with a higher Y. There are GC with multiple 
subgiant bran ches (e . g. M2 2 (jMilone et al.ll2010al lbh. NGC 
104, etc.). See lPiottd (|2009l ') for a review of the evidence. 

These second generation stars form a substantial por- 
tion of all stars in GC ( up to 60-70 percent of the total 
number of stars in a GC (|Carretta et al.ir2009l )l and should 
have a substantial impact on what is observed during the 
AGB and PN phases of evolution. In this paper I look at 



2 MODELS 

Most of the relevant det ails of this model are ex- 
plained in iBuell et all (Il997l ). iBuelll l| 19971 ). and 
iGavilan. Buell. fc Mollar(|2005l '). In this section I con- 
centrate on modifications of this model. Of particular 
importance is the mass-loss on the RGB and E-AGB which 
can be significant input, especially at the low ZAMS masses 
of globular cluster stars. Particular attention is paid to the 
effect of enhanced helium abundances on mass-loss during 
these stages. The mass-loss in these stages is found by in- 
tegrating the mass-loss rate formula over the Padua tracks. 
There is an important point to note, in the pre-TP-AGB 
mass-loss model the stellar evolution and the mass-loss are 
not coupled and hence the equations for mass-loss derived 
below should be used with care. The mass-loss shift is 
probably a small effect as explained below but the reader 
should be aware of it. 



2.1 Red giant mass-loss 

The mass-loss which occurs on the red giant branch (RGB) is 
very important for low mass stars found in globular clusters, 
in some extreme cases it may sometimes prevent the star 
from even reaching the TP-AGB. Most of the pre-TP-AGB 
mass-loss occurs during the RGB. The standard method to 
determine the amount of mass-loss is to use Reimers' Law 
(|Reimerslll975h given by 

M = ,— (1) 

where L, R and M are the stellar luminosity, radius and 
mass, respectively in solar units. However, to calculate the 
pre-TP-A GB mass-loss the modified version of the Reimers 
formula of ISchroder and CuntjIioOSi given by: 



M MOOOK' 



'(1 + 



P0 



(2) 



4300p* ' 

where Tofi is the effective stellar temperature, is the sur- 
face gravity of the star in cgs units. Values of 27400cms~^ 
for (jQ and 8.0 x 10~ for r) were ad opted (Value recom- 
mended by I Schroder and Cuntzll2005h . This new mass-loss 
rule appears to give better re sults for horizonatal bran ch 
masses then the Reimer's rate (|Schroder and Cuntz|[2005l '). 

This mass-loss law is applied to the variable Y stel- 
lar evolution tracks from the Padova stellar evolution- 
ary library (http://plcadi.pd.astr o.jt| described in detail in 
iBertelli et al.. (.20081 and .Bertelli et all (|2009l '). To determine 
the red giant mass-loss the mass-loss rate was integrated 
from the beginning of the red giant branch (encoded in the 
Padova files as brgbs) up to the tip of the red giant branch 
(encoded as trgb) using the trapezoidal rule. The amount of 
mass-loss between time steps in the models is given by: 



AA/i = -(fm. 



-\ + rh,)(t 
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(3) 
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where ti and ti+i are the model times and rhi+i and rhi are 
the mass-loss rates at the corresponding times. The total 
mass-loss is determined by summing all of the AMiS. 

Table [T] shows the values of Y and Z for which mass- 
loss was computed. For all available masses the mass- loss on 
both the EAGB and RGB was computed. Z is the value of 
Z on the ZAMS and Y is the value of Y on the ZAMS. 

The mass-loss on the RGB as a function of the ZAMS 
mass for all available values of Y and Z is shown in figure [1] 
In all panels it is evident the amount of mass-loss decreases 
as Y increases. This occurs because stars with higher values 
of Y means the RGB star will have smaller radii and higher 
surface gravity due to the lower opacity in the outer layers. 
These factors lower the mass-loss rates and the total mass- 
loss. 

The RGB mass-loss was fit using two linear fits for 
higher and lower ZAMS masses. The transition point be- 
tween the fits was determined by visually estimating the 
mass where the slope appears to change. This mass is typi- 
cally found around a ZAMS mass of 0.8-0.9 Mq. The higher 
mass fit was terminated where the high mass line crosses the 
horizontal axis. This termination point was also estimated 
visually. For masses larger the mass where the higher mass 
fit crosses the axis the mass-loss is 0. The equations of the 
linear fits for low and high masses are given by AMrgb,1ow 
and AMRGB.high. To make the fits work at the two highest 
metallicities {Z — 0.008 and 0.017) the fitting was done by 
excluding the M = 0.6 Mq models. These can be safely ex- 
cluded since the mass-loss for these very low mass stars elim- 
inates their envelope before the tip of the RGB is reached 
and such models will not be considered in this paper. Visual 
inspection indicates the fits are in good agreement to the 
mass-loss calculations. 

The equations for AMrgb,1ow and AMRGB,high are 
given by 

AMrgb.Iow ^aiiM + aio (4) 

AMRGB,high = a2iM + a2o (5) 

The mass-loss is found by calculating the value of both AMs 
and finding the maximum value. If the mass-loss is found to 
be negative then the value of the mass-loss is set to 0. The 
coefficients of these equations for the different values of Y 
and Z are shown in table [2] 

No attempt has been made yet to calibrate this mass- 
loss, which will be done in a later paper. However, the mass- 
loss values from these equations appear to be reasonable. For 
example a 1.0 Mq Y = 0.26 Z — 0.017 star would experience 
0.28 Mq of mass-loss on the RGB which is typical of other 
models. A typical turn-off mass of 0.80 Mq with Y = 0.245 
and Z = 0.0008 gives a RGB mass-loss of 0.22 Mq which 
is reeisonable giving a zero-age horizontal branch mass of 
appr oximately 0.58 MQwhich is similar to measured values 
(e.g. iGratton et aD (j2010l ll. 

It should be noted, as suggested by the referee, that the 
method used to find the mass-loss is not consistent with the 
stellar evolution models. As the star loses mass its surface 
gravity would decrease causing the star to expand. For the 
models used this would result in a higher mass-loss rate near 
the tip of the RGB and a greater amount of mass-loss on 
the RGB (and the E-AGB) then is calculated here. However, 
this effect should be relatively small since the deviation will 



only be really significant at the tip of the RGB. Although 
the method used here is not strictly consistent the relative 
differences in mass-loss due to the effect of the ZAMS helium 
abundances and the ZAMS metallicity should be correct. 



2.2 First dredge-up 

This is important since when all stars enter the RGB the 
convective envelope penetrates into regions of the star where 
partial H burning has occurred bringing these products up 
to the surface. If a PN is formed this process will have mod- 
ified the surface abundances. This paper only considers the 
effects on He and the CNO elements since these are what 
is observed in planetary nebula. The first dred ge-up (FDU) 
prescription of iGroeneweeen fc de Jon3 (j 19941 ) is used. 



2.3 Mass- loss on the early- AGB 

The same procedure was applied to the early- AGB (E-AGB) 
portions of the Padova tracks. An additional condition of 
starting the mass-loss when the temperature was below 
4500K was assumed since this mass-loss law is applicable 
only to K and M stars. 

Figure [2] shows the calculated mass-loss during the E- 
AGB and the fits to these mass-losses. The most obvious 
trend is as Y increases so does the amount of mass-loss. 
This is opposite to the trend on the RGB. In this case the 
core mass on the AGB is increased which also increases the 
luminosity. This increase in the luminosity on the E-AGB 
results in greater mass-loss. This enhancement of mass-loss 
on the E-AGB is important since it means the higher the 
value of Y the more mass is lost on the E-AGB. It means 
such a star has a higher probability to reach the horizontal 
branch but its envelope may not survive to reach the first 
thermal pulse. 

The mass-loss on the E-AGB is fit using 4 fits in differ- 
ent regions of mass. The lowest mass range (M < 1.5 Mq) 
is fit via a cubic, the next mass range up (1.5 Mq < M < 
2.0 Mq) is fit using a quadratic fit. The next mass range up 
(2.0 Mq < M < 4.5 Mq) is fit using a linear fit. Finally the 
highest masses are fit using a constant value of mass-loss. 
The points of intersection between adjacent fits were visu- 
ally estimated. This procedure gives a good fit to the model 
mass-losses. 

The equations for the E-AGB mass-loss in the first two 
mass regions are given by: 

AMe-agb = aislVf + ai2M^ + auM + aio (6) 
and 

AMe-agb = a22M^ + a2iM + 020 (7) 

where M is the mass of the star on the ZAMS. Only the co- 
efficients of first two regions have been included in table|3]to 
save space and since no models of sufficient mass which need 
the fits for the upper regions are calculated in this paper. 
The E-AGB mass-loss is calculated by finding the intersec- 
tion of the two regions and then choosing the appropriate 
region and plugging into the corresponding equation. 
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Table 1. Padova Models Used 

Z Y 



0.0001 0.23,0.26,0.30,0.40 

0.0004 0.23,0.26,0.30,0.40 

0.001 0.23,0.26,0.30,0.40 

0.004 0.23,0.26,0.30,0.40 

0.008 0.23,0.26,0.30,0.40 

0.017 0.23,0.26,0.30,0.34,0.40 
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Figure 1. Each panel in the figure shows the calculated mass-loss on the red giant branch for the Z=0.0001, 0.0004, 0.001, 0.004, 0.008 
and 0.017 models for different available values of Y. The open squares, open circles, open triangles and closed squares are the mass-losses 
for y=0.23, 0.26, 0.30 and 0.40, respectively. The upside down triangles in the Z = 0.017 panel are the mass-losses for the Y = 0.34 
models. The solid, dashed, dotted, and dot long dashed lines are the fits for the y=0.23, 0.26, 0.30, and 0.40 mass-losses, respectively. 
The dot dashed line is the fit for the y=0.34 models in the Z=0.017 panel. 



2.4 Core mass at the first pulse 

From the Padova stellar evolution library I also extracted 
the core mass as a function of the mass. Figure [3] shows the 
mass of the carbon-oxygen core as a function of mass for 
several values of Y and Z. An important point to note is 
that as the initial helium mass fraction increases so does the 
mass of the core. The core mass is important since it is the 
most important factor controlling the luminosity of an AGB 
star. 

Figure [3] show the core mass at the first pulse from the 



Padova models for Z =0.0001, 0.0004, 0.001, 0.004, 0.008, 
and 0.017 and Y =0.23, 0.26, 0.30 and 0.40 as a function of 
mass. Each set of models with a given Y and Z have been fit 
by a double quadratic fit, one at lower masses < 1.5 Mq and 
one at the higher masses. The transition between the two 
was found between 1.3 and 2.0 Mq. The transition between 
the lower mass and higher masses was determined by visual 
inspection of where the core mass begins to rise steeply. In 
all cases the fits to the points are good giving core masses 
which are typically less than 0.02Mq difference. 

The equations of the quadratic fits for low and high 
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Table 2. Red giant mass-loss coefBcients 



Y 




ail 








0121 


0120 


0.0001 


0.23 


-0.53375 


0, 


.587959 


-0.146583 


0.257894 


0.0001 


0.26 


-0.48794 


0, 


.528645 


-0.158946 


0.263457 


0.0001 


0.3 


-0.461671 


0, 


.509082 


-0.128568 


0.210224 


0.0001 


0.4 


-0.317215 


0, 


.330223 


-0.0682917 


0.0996025 


0.0004 


0.23 


-0.556562 


0, 


.651869 


-0, 


.171144 


0.304073 


0.0004 


0.26 


-0.54755 


0, 


.618164 


-0, 


.173845 


0.298794 


0.0004 


0.3 


-0.50535 


0, 


.560052 


-0, 


.166315 


0.275781 


0.0004 


0.4 


-0.52254 


0, 


.499175 


-0, 


.142323 


0.210173 


0.001 


0.23 


-0.607168 


0, 


.728346 


-0, 


.197265 


0.353301 


0.001 


0.26 


-0.571589 


0, 


.678889 


-0, 


.161241 


0.298386 


0.001 


0.3 


-0.509972 


0, 


.604365 


-0, 


.190175 


0.317975 


0.001 


0.4 


-0.384763 


0, 


.443838 


-0, 


.170916 


0.258054 


0.004 


0.23 


-0.56752 


0, 


.764646 


-0, 


.223724 


0.418548 


0.004 


0.26 


-0.65305 


0, 


.815914 


-0, 


.234945 


0.419982 


0.004 


0.3 


-0.658095 


0, 


.794849 


-0, 


.236612 


0.403639 


0.004 


0.4 


-0.45996 


0, 


.550102 


-0, 


.234695 


0.359844 


0.008 


0.23 


-0.73907 


0, 


.973301 


-0, 


.252334 


0.479442 


0.008 


0.26 


-0.75636 


0, 


.960051 


-0, 


.263265 


0.478337 


0.008 


0.3 


-0.831625 


0, 


.984962 


-0, 


.259694 


0.453806 


0.008 


0.4 


-0.629659 


0, 


.728283 


-0, 


.262052 


0.40875 


0.017 


0.23 


-0.673745 


0, 


.971984 


-0, 


.255924 


0.50791 


0.017 


0.26 


-0.682876 


0, 


.959522 


-0, 


.241278 


0.465487 


0.017 


0.3 


-0.75486 


0, 


.982077 


-0.27485 


0.495803 


0.017 


0.4 


-0.7626 


0, 


.887937 


-0.280617 


0.448454 


0.017 


0.34 


-0.6274 


0, 


.826845 


-0.256405 


0.452657 



Table 3. Coefficients for fits to early-AGB mass-loss 



z 


Y 


ai3 


ai2 


ail 


0.0001 


0.23 


-0.0217193 


0.0915006 


-0.122712 


0.0001 


0.26 


-0.215429 


0.683192 


-0.699251 


0.0001 


0.3 


-0.00119991 


-0.0408501 


0.10212 


0.0001 


0.4 


-0.0758491 


0.117157 


0.0797828 


0.0004 


0.23 


-0.0938861 


0.306061 


-0.327053 


0.0004 


0.26 


-0.0451582 


0.138365 


-0.144856 


0.0004 


0.3 


-0.0753496 


0.168929 


-0.0967276 


0.0004 


0.4 


0.174151 


-0.704201 


0.939363 


0.001 


0.23 


-0.0212184 


0.0804842 


-0.116847 


0.001 


0.26 


-0.0886186 


0.296724 


-0.336025 


0.001 


0.3 


-0.129515 


0.334909 


-0.262369 


0.001 


0.4 


0.288478 


-1.1637 


1.43381 


0.004 


0.23 


0.00637805 


-0.00833866 


-0.0338022 


0.004 


0.26 


0.00369947 


-0.00247321 


-0.0397723 


0.004 


0.3 


-0.0537507 


0.173439 


-0.212193 


0.004 


0.4 


-0.689675 


1.65297 


-1.08412 


0.008 


0.23 


-0.0179112 


0.0872929 


-0.148306 


0.008 


0.26 


-0.00338358 


0.0399171 


-0.111972 


0.008 


0.3 


-0.0206433 


0.0880253 


-0.162452 


0.008 


0.4 


1.46767 


-5.9136 


7.68436 


0.008 


0.34 


-0.012937 


-0.00120956 


-0.0132868 


0.017 


0.23 


-0.00576822 


0.041759 


-0.106142 


0.017 


0.26 


-0.053432 


0.248533 


-0.399333 


0.017 


0.3 


0.089118 


-0.258045 


0.170682 


0.017 


0.4 


-0.149333 


0.461895 


-0.471668 


0.017 


0.34 


0.399763 


-1.37573 


1.41367 



masses are given by Mco,iow and Mco.high- The equations 
are: 

M cO.low — 0tl2 + aiiM -I- aio (8) 
Mco.high = a22M^ + 021 M + a2o (9) 
where M is the ZAMS mass. The coefRcients for the different 



aio 


a22 


a2i 


a20 


0.0744702 


0.130361 


-0.514963 


0.508129 


0.256586 


0.0581723 


-0.180144 


0.136288 
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values of Y and Z are shown in tabled The procedure used 
is to find the point of intersection between the two fits and 
then to plug in the relevant mass. 

The most obvious trend is there is an increase in the core 
mass as the value of Y increases. This is important since on 
the AGB a larger core mass leads to a higher luminosity. 
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Figure 2. Each panel in the figure shows the calculated mass-loss on the early asymptotic giant branch as a function of ZAMS mass 
for the Z =0.0001, 0.0004, 0.001, 0.004, 0.008 and 0.017 models for difl'erent available values of Y. The symbols and Unes have the same 
meaning as those in figure 1. 



This is also important since the mass at the first pulse is an 
important factor in determining the mass of the CSPN. 



2.5 Third dredge up 

In synthetic AGB models the standard method to model the 
third dredge up (TDU) effect is to use a dredge-up param- 
eter A so that 



A 



dredge 



(10) 



where AMdrodgo is the mass dredged up and AMc is the 
increase in the core mass during the preceding interpulse 
phase. During a thermal pulse the star develops a convective 
shell in the region between the intershell region between the 
base of the hydrogen-rich envelope and just above the core. 
This region is helium and carbon rich since it consists of 
the products of partial helium burning. At the end of the 
thermal pulse the convective envelope may penetrate into 
this region and mix this carbon and helium rich material into 
the envelope. The parameter A is a measure of how deeply 
the convective envelope penetrates into this intershell region 



and determines how much mass is mixed up into the outer 
layers. 

A number of authors have used synthetic models to con- 
strain the value of A and the mini mum core mass value at 
which TDU can occur, Mc.min (e.g. iGroenewegen fc de Jon^ 
(|l993MMarigo. Girardi. fc Bressanll 19991 )). These authors 
used the LMG and SMC carbon star luminosity functions 
to constrain both Mc.min and A parameters. They show the 
value of Mc.min is approximately 0.58 Mq. In this paper the 
value of Mc,miii is treated as a free parameter but values near 
this are always chosen. The value of A is a free parameter 
but it will always be small since most models will experi- 
ence only one TDU and I expect the value of A would grow 
during subsequent thermal pulses if they were to occur. 



2.6 Mass-loss on the TP- AGB 

On the TP-AGB, ma ss-loss is calculated by the p ulsation 
period-mass loss law of lVassiliadis and Wood ((1993) without 
their correction for periods above 500 days. To make the 
transition from the modified Reimer's rate to this pulsation 
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Figure 3. The figure shows the core mass at the onset of the first pulse for the Z =0.0001, 0.0004, 0.001, 0.004, 0.008, and 0.017 models. 
The symbols have the same meaning as they do in figure 1. 



Table 4. Coefficients for the fits to the first pulse core masses 
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mass-loss rule the modified Reimer's rate is used until the 
pulsation mass- loss rule becomes larger. 



2.7 Conditions for formation of a visible 
planetary nebula 

The formation of a visible PN requires that the central star 
of the planetary nebula (CSPN) reach a temperature of ap- 
proximately 30000K on its blueward journey from the AGB 
to the white dwarf cooling tracks before the ejected envelope 
has sufficient time to disperse into the interstellar medium. 
In a globular cluster the typical star at the turn-off would 
produce a CSPN of mass approximately 0.52 Mq. 

The criteria used in this paper assumes that if a plan- 
etary nebula does not form before the ejected envelope ex- 
pands to 0.5 pc it will not be visible. If an expansion rate of 
15 km • s~^ is assumed for PN then the maximum transition 
time to form a visible planetary nebula is given by approxi- 
mately 25000 years. Since PN with expansion ages of 30000 
years exist this conservative value is adopted. 

To find the approximate transition time for all the Z — 
0.016 models of IVassiliadis and WoodI l|l994 ) were linearly 
interpolated in time and logTefi to find the time when the 
effective temperature becomes 30000K. This is the transition 
time ttrans- lu figurc [J] logttrans is plotted as a function of 
the CSPN mass. The results for the Z = 0.016 models are 
fitted with both a quadratic and cubic fit. These fits are used 
to extrapolate an approximate range of values for CSPN 
masses lower than 0.57 Mq cores. Also included were the 
results from other metallicities. With the exception of one 
point from the Z = 0.004 models (their M = 2.0 Mq model) 
all of the points follow the trend of the cubic. 

The quadratic fit is given by: 

logitrans = 19.3316M^^ - 36.0874Mc 18.0464 (11) 
The cubic fit is given by: 

/opttrans = - 181 .893Afc -^41 7.538Mc - 322.289A/c +85.6( 12) 

The cubic fit is a better fit in that it passes closer to all 
the points than does the quadratic fit. The cubic fit seems 
to better capture the extrapolation, however, since these 
are used for extrapolation caution needs to be exercised 
and the values derived here should be treated as approxi- 
mate. If the transition time is set to 30000 years and the 
quadratic and cubic equations are solved for Mc the results 
are O.522M0 and 0.548Mq, respectively. The cubic extrap- 
olation is weighted more heavily and the adopted minimum 
value of the CSPN mass to produce a visible PN is 0.545Mq. 
This value is very close to the lowest observed value of a 
globular cluster CSPN central masses, JaFu 1 and JaFu 2. 
It should be noted a choosing a different set of post- AGB 
models might lead to different conclusion about the mini- 
mum CSPN mass because the transition time in post- AGB 
models depends strongly on the adapted mass-loss rate. 



2.8 TP-AGB models 

The TP-AGB is followed using a synthet ic AGB code which 
is a descendent of the iRenzini fc Volil (119811 ) code. The 



code begins with a guess at Tefi to calculate the surface 
boundary condition. The equations of stellar structure is 
then integrated to the base of the convective envelope. The 
value of the effective temperature is modified until the base 
of the convective envelope is at the same position as the 
core mas s. The opac i ties u sed for high temperatures are 
the laicsi as fc Rogers! lll996l) and for low temperatures the 
opacities of Alexander fc FergusonI (|l994l ) are used. The lu- 
minosity of the star is calculated u sing the expressions in 
IWagenhuber fc GroenewegenI l|l998l ). A mixing length pa- 
rameter, a = i/Hp, of 1.70 is used. This value is chosen 
since it is close to typical values of values of a chosen for 
solar models. 



3 RESULTS 

3.1 First Generation Star Models 

First generation stars are stars with the lowest possible Y of 
a globular cluster. A series of models with masses between 
0.7 and 1.0 Mq were calculated for metallicity values of 
0.0002 and 0.004 and a typical Y value of 0.245. Figure [5] 
shows the final white dwarf mass as a function of the ZAMS 
mass for a variety of Z values and the Y value of 0.23. If 
we adopt a typical globular cluster age of ~ 13 Gyr it gives 
typical turn-off masses of ~0.85 Mq. Note for all values of 
the metallicity, Z, the final mass is below the mass needed 
to produce a visible PN. 

If these stars produced visible PNe then there would 
be a number of PNe in globular clusters with He/H around 
0.100 and values of N/0 and C/O similar to the values of a 
typical disk PN. Since the CSPN masses are too low, these 
models suggest first generation stars do not produce visible 
PNe as the ejected mass would disperse before it could oc- 
cur. This reduces the number of expected PNe in GCs since 
these stars are too low mass to produce them. The larger 
the fraction of primordial stars is the lower the number of 
expected PNe. Since no typical PN are observed in the GC 
system this model explains this observation. 

The typical white dwarf mass which would result from 
these first generation stars is about 0.525 Mq. These values 
are in rough agreement with the typical values inferred for 
the WD mass at the top of the coo l ing se quence determined 
from measurements. [Hansen et al.l l|2007t ) finds WD masses 
at the top of the WD cooling sequence of globular cluster 
NGC 6397 between 0.50 an d 0.53 M^ although their error 
analysis favors a lower value. [Hansen et al.l (|2004l ) found the 
mass of the WD at the top o f the cooling sequenc e of glob- 
ular cluster M4 as 0.55 Mq. IStrickler etal] (120091 ') find the 
mass of the WD cooling sequence of CO WDs in NGC 6397 
is 0.53 Mq. 



3.2 Planetary Nebulae from Second Generation 
Stars 

These stars start on the ZAMS with a high initial He abun- 
dance (y ~ 0.30). Figure [6] shows the calculated initial-final 
relationships for ZAMS stars with masses ranging from 0.7 
to 1.0 Mq with Z = 0.004 and a range of possible values 
of Y . It shows as the initial helium abundance is increased 
the final white dwarf mass increases. This effect is somewhat 
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Figure 4. The figure show s the transition times as a function of core mass. The open squares are the results calculated from the 
IVassiliadis and WoodI l ll994) Z = 0.016 models. The filled squares, open circles and filled circles are the results from the Z =0.008, 0.004 
and 0.001 models, respectively. The solid and dashed lines are the quadratic and cubic fits described in the text. 
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Figure 5. This figure shows the model CSPN masses as a function of ZAMS mass. The open squares are the Z = 0.0002 models and 
the closed squares are the Z = 0.004 models. 



counter balanced by the shorter ZAMS lifetimes of stars and 
lower masses of the turn-off with higher initial Y . However 
even for an assumed turn-off of 0.75 Mq with the higher val- 
ues of Y the final mass is near or above the limit to produce 
a visible PN. 

For the stars with a primordial value of helium (Y = 



0.245), Z = 0.004 and a turn-off mass of 0.85 Mq the core 
mass here is 0.529 Mq which would be too small to produce 
a visible PN. These models suggest the reason only PN with 
high He/H appear is this allows larger core masses; which 
lead to shorter transition times allowing a visible PN to 
appear. 
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Since these He enhanced stars can be a significant frac- 
tion of the total number of stars in a globular cluster it also 
predicts that there where two separate populations with dif- 
ferent values of Y exist there should should be two CO WD 
cooling tracks at slightly different masses. Since it appears 
the typical WD mass is around 0.53 Mq and 0.54-0.55 Mq 
we estimate that these tracks are separated by 0.01-0.02 Mq. 
This would be observationally challenging to do but might 
be possible. 

Figure [7] shows the model grids used to fit the CSPN 
masses and He/H for the model parameters which are the 
closest fit to observed values of A/cspn, He/H and log O/H-f 
12 from Jacoby et al. of JaFu 1 and JaFu 2. For both nebula 
the ZAMS values of Z was varied until the model value of 
0/H was a close match to the observed abundances of 0/H 
found in the nebula (from Jacoby et al.). Then the value of 
Y was adjusted until a close fit to He/H was obtained. The 
goal here was to find a reasonable set of parameters and not 
to fine tune this to the best possible fit. 

For JaFu 2 the best fitting model has a core mass of 0.55 
Mq which is the same as the observed value. This core mass 
is produced by a model ZAMS star of mass approximately 
0.85 Mq which is a quite reasonable turn-ofi' mass. The val- 
ues of the model ZAMS Y and Z are 0.282 and 0.00102, 
respectively. The value of Y is quite reasonable in this case 
being in line with which are quite reasonably in line with 
expectations for helium enhancements. JaFu 2 is found in 
globula. r cluster NGC 6441 which has an [Fe/H] of -0.46 
l|Harrisl (|l996 - 2010 editioni )), which is a bit higher than 
the adopted value of Z for the nebula. However, [Fe/H] is 
not actually fit here and O/H is. It is quite possible there are 
oxygen depletions which correspond to the helium enhance- 
ment here over the different populations in this cluste r . This 
would be similar to what was found bv lPiotto et~al] ()2007l l 
for NGC 2808. It is expected He is enhanced by the opera- 
tion of the CNO cycle which would result in a depletion of 
oxygen. 

For JaFu 1 the best fitting model has a final CSPN mass 
of 0.54 Mq, Fzams = 0.333, Zzams = 0.006947. The value 
of Y is typical of the expected Y enhancements. However 
the ZAMS mass of the "best" model is 0.70 Mq which is 
too small but the errors are large enough here to allow a 
range of possible ZAMS masses up to 0.90 Mq. The ZAMS 
mass of 0.90 Mq gives a core mass of 0.558 Mi^ which is 
a littl e bigger than the observed value from Ijacobv et al.l 
l| 19971 ') (0.54 ± 0.2 Mq) but it is well within the error bars. 

It is difficult to compare the metallicity of the model for 
JaFu 1 to that of the cluster Pal 6. A range of potential val- 
u es are found in the lite ratur e. On the lo w end are the works 
of iLee fc Car"nevl (|2002l ) and lZinnI (|l985l ) who find [Fe/H] of 
— 1.08±0 .06 and -0.74, resp ectively. On the high end are the 
works of 'Bica ct al.' (1998) who find [Z/ZqJ = -009 and 
lOrta lani. Bica & Barbuy (1995) and iMinnitil (j 19951 ') who 
find [Fe/H]=-0.4 and -1-0.22, respectively. The oxygen abun- 
dance of JaFu 1 suggests an overall metallicity between 1/3 
and 1/2 solar which is within the metallicity range on the 
high end of the metallicity possibilities. 

The observed nebular mass of JaFu 2 (of 0.04 Mq from 
Jacoby et al.) is a reasonable match to the models. Low 
ZAMS mass models with values of Y around 0.30 predict 
a nebular mass is just a few times 0.01 Mq. However, the 
observed nebular mass of JaFu 1 is much too high (0.40 Mq) 



but Jacoby et al. determined this using an assumed filling 
factor of 1 and note a filling factor of 0.1 gives a much lower 
mass. 

These models provide a good fits to the most important 
parameters; the CSPN mass, the observed value of He/H, 
the turn-off masses and the oxygen abundance of these PNe. 
These two nebula can be explained with reasonable progen- 
itors and slightly enhanced values of He/H. 

3.3 High C/O Nebula 

No star with a mass below the turn-off of any globular clus- 
ter will produce a carbon-rich planetary nebula like Ps 1 
since a core mass of approximately 0.58 Mq is required for 
the third dredge-up to occur and since globular cluster stars 
at the turn-off will have core masses well below this. There 
is however a channel to get to these masses - blue stragglers 
are thought to be merged stars and as such are more massive 
than stars on the turn-off. Since blue stragglers in principle 
could have masses up to twice that of the turn-off these are 
the most likely channel. 

The model u sed here is the same as discussed in 
iBuell et all (|l997l ) where after a dredge-up occurs the en- 
velope is so polluted with ca rbon it is es sentially immedi- 
ately ejected. As discussed bv lBuell et al ] ([1997) this single 
dredge-up heavily pollutes a low metallicity envelope and 
the metallicity of the envelope jumps to near solar values. 
This results in an expansion of the star and this resulting 
expansion causes the star to immediately switch into a su- 
perwind phase. The remaining envelope is eliminated in a 
few thousand years suggesting the carbon star phase would 
be shorter than the PN phase. This explains the lack of car- 
bon stars in globular clusters. 

The models calculated here for low metallicity (Z ~ 
0.0002) confirm this basic scenario. If Mc,min is set to 
0.57Mq then the minimum ZAMS mass is 1.15 Mq. Fig- 
ure[8]shows the TP-AGB evolution of the important param- 
eters of this model. On the last pulse the envelope is heavily 
polluted and note the resulting change in the mass-loss rate. 
The envelope in this case is ejected in a few thousand years. 
The model produce a CSPN of the correct mass and the 
C/O ratio of the model compares favorably to observations 
of the final ejecta. 

Since Ml 5 co ntains a number of blue stragglers 
ijYannv et al.l (| 19941 )'). which, in this model are the progeni- 
tor star for a high C/O nebular like Ps 1, M15 is a reasonable 
site for the formation of this PN. 

The hardest observation to match is the low nebular 
mass of Ps 1. In the blue straggler model the predicted neb- 
ular mass is ^ 0.2 Mq whereas the observed nebular mass is 
about 1/lOth this value. The observed smaller nebular mass 
suggests for the blue straggler model to be a viable model 
their must have been additional mass-loss, perhaps due to 
binary interaction. There seems to be no way to reproduce 
this observation with a single star. 

An alternative scenario is described by lOtsuka et al.l 

who argue stars like K648 have evolved from a bi- 
nary and its progenitor is a CEMP-s star. The problem for 
this scenario is to match the CSPN mass of K648 the ZAMS 
mass of the star needs to be larger then the turn-off mass. 
Otherwise the CSPN mass would be smaller. 

Another complication to the model scenario in this pa- 
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Figure 6. Figure shows the model CSPN masses as a function of ZAMS mass for models with Z = 0.004. The open squares, filled 
squares, open circles, and filled circles are the Y =0.244, 0.274, 0.304, and 0.334 models, respectively. 
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Figure 7. This figure shows H!c/tl as a function of CSPN mass for JaFu 1 and JaFu 2 as filled squares. The models of 'best' fit for JaFu 
1 and JaFu 2 discussed in the text arc shown as open circles and open squares, respectively. Each of the points on the model curves 
starts at 0.70 Mq and is iterated by 0.05 Mq up to 1.00 Mq. 



per is if the initial metallicity of the star is increased the 
minimum mass to produce a carbon star gets larger. This 
occurs because the mass of the core at the first pulse gets 
smaller with smaller metallicity. This is consistent with stud- 
ies which infer the minimum mass to produce a carbon star 
in the SMC, LMC and the Galaxy where as the metall- 



icy goes up so does the minium mass (See - references here 
(Marigo, Groenewegen, etc.)). In this model a blue straggler 
in a higher metallicity cluster {Z = 0.004) would be less 
likely to produce a carbon-rich nebula. This blue straggler 
would produce a visible PN with an oxygen-rich composi- 
tion. 
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Figure 8. This figure shows the TP-AGB evolution of a ZAMS star M=1.15M0, Z = 0.0002 and Y = 0.2454. The top panel shows the 
evolution of the ratio of carbon to oxygen. The lower panel shows the evolution of the mass-loss rate (MLR). The points are the values 
at the end of each interpulse cycle. 



Blue stragglers are rare but prominent members of GCs 
and cannot be expected to produce many PNe. A simple es- 
timate can be obtained in the following manner; assume the 
average lifetime of a blue straggler is 10® years and the life- 
time of a PN is 25000 years. If it is assumed every blue strag- 
gler produces a PN then we should expect 1 PN for every 
40000 blue stragglers. If we assume the number of blue strag- 
glers per cluster is 100-200 and use 150 for the number of 
globular clusters then the total number of blue stragglers in 
the GC system should be between 15000-30000. This would 
give the number of PN from this channel as 0.375-0.70 PNe. 
The assumed numbers are probably optimistic since the life- 
time of blue stragglers could be larger; the lifetime of PN 
smaller and the 100-200 is the number of blue stragglers for 
bigger clusters. The actual number of expected PN may be 
even smaller. If it is assumed there are sufficient numbers to 
produce only one in the entire GC system then it is possible 
Ps 1 is the only example and an oxygen-rich PN is equally 
likely. 

This model shows it is possible to produce carbon-rich 
PNe in globular clusters but it is not possible for it to match 
the nebular mass of Ps 1 , which can only be done by invoking 
an ad hoc binary explanation. This may be reasonable since 
this star probably formed in a binary but at the moment it is 
untestable since there are no similar PN in the GC system to 
compare it to. A similar PN with a well constrained CSPN 
mass and nebular mass would tell us if K648 is typical of its 
class. The only possible known analogs are the halo PN such 
as BB-1 and H4-1 which have C/0>1 (see Howard et al.). 
Unfortunately the distances to these nebula is poorly known 
and it will be difficult to get precise information about the 
CSPN and nebular masses for comparison purposes. 



3.4 What is the origin of the planetary nebula in 
M22? 

M22 contains a ver y interesting PN GJJC-1 which is 
nearly hydrogen free l|Gillett et al.|[l989h because no neb- 
ular hydrogen lines have been detected. The nebula ap- 
pears to be overabundant in O and Ne. The central star 
of G JJC-1 has been observed (|Harrington and Paltogloul 
Il993l ) and it is suggested that the surface abundances are 
He/H=0.5, C is 6 times solar a nd N is 14 times solar 
l|Harrington and Paltogloul (| 19931 ) . iRauch et al.l (119981 )). 

There are only two other hydrogen free nebula out of 
approximately 1000-2000 total known in the Galaxy. Given 
the inferred percentage of a PN being hydrogen free of 0.1- 
0.2 percent it seems unlikely that this PN simply occurred by 
chance. One speculative possibility is the existence in some 
clusters of a third generation of stars which would have very 
extreme {Y = 0.40) helium enhancement. 

Recent work on th e cluster by iDa Costa et al.l (|2009| ) 
and lMarino et ahl (|2009| ) indicates M22 has at least two sep- 
arate populations which vary in their value of [Fe/H]. The 
average [Fe/H] is between -1.7 and -1.9 but it ranges up to 
-1.4. The higher [Fe/H] population probably has a higher 
value of Y . It has been suggested M22 is the remains of 
a captured dwarf elliptical galaxy. As such it seems a rea- 
sonable candidate for a population with extreme helium en- 
hancement. M22 may also have a more complicated star for- 
mation history with the formation of different populations 
spaced out by significant amounts time. 

A model with M = 0.80 M©, Y = 0.40865, and Z = 
0.000886 was run. In this model the star reached the hori- 
zontal branch but the entire envelope was lost by the end of 
the E-AGB. The mass-loss on the E-AGB increases as the 
value of Y goes up because the core mass goes up. The lumi- 
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nosity of an AGB star will depend strongly on the core mass 
and this increased luminosity drives the increased mass-loss. 
This model produces a core mass of 0.60 Mq. If the mass is 
lowered to 0.75 Mq then the model once again produces a 
HB star but all the remaining envelope mass is lost on the 
E-AGB. This lower mass produces a 0.58 Mq remnant. 

Why is this significant? This is close to the inferred 
mass of GJJC-l's CSPN and production of a hydrogen free 
nebula requires getting rid of the initial hydrogen before the 
PN phase. Getting rid of all the envelope means if this star 
experiences a thermal pulse then there would be little to no 
hydrogen ejected and the star would eject helium rich ma- 
terial and it would look much like GJJC-1. This suggestion 
should be regarded as speculative and the positive evidence 
for it is thin but the idea seems to be possible and merits 
additional study. 

3.5 How many of each type? 

Without a detailed population study which is beyond the 
scope of this paper it is only possible to show the numbers 
work out approxim ately correct. I start by assuming the 
iJacobv et aLfl|l99'iD statement that given the total luminos- 
ity of the GC system then the number of PNe should be 16. 
If it is assumed that 70 percent of globular cluster stars are 
primordial and have a turn-off in the 0.80-0.90 Mq range 
then this reduces the number to about 4-5 since these stars 
produce none. The number from this appears to be 3 which 
is in rough agreement. Further, assuming all blue stragglers 
produce a visible PN, we assume about 1 from this part 
which matches the designation of Ps 1 (K648) as being pro- 
duced by this channel. The rest will be produced by second 
generation stars (which may or may not all produce a visible 
PN) . Therefore the numbers are roughly consistent with the 
number of known PNe. 



4 DISCUSSION 

What do JaFu 1 and JaFu 2 tell us about second generation 
GC stars? They confirm, independently of the use of colour 
magnitude diagrams, stars exist in globular clusters with 
Y = 0.28 — 0.33 and it suggests these second generation 
stars are a significant fraction of the number of GC stars. 
These two PNe could turn out to be very important since 
they allow the direct observation of elements which can not 
be observed directly using stellar spectroscopy. 

For JaFu 1 log N/0 = -0.52 (Jacoby et al.) 
which is slightly higher t han logN/Og — —0.88 
l|Asplund. Grevesse. fc Sauvall |2005| ) . This ratio is consis- 
tent with both a first dredge-up event and possibly a small 
amount of nitrogen enrichment or oxygen depletion. For 
JaFu 1 logS/0 = —1.35 which is consistent with the so- 
lar logS/O — —1.50. From these abundance ratios either 
all these have been enriched by the same relative amount 
or the star has not been enriched relative to a first genera- 
tion star in Pal 6. The sulfer abundance of JaFu 1 is given 
by logS/H — logS/HQ — —0.55. This is consistent with a 
cluster with metallicity between 1/3 and 1/2 solar. If the 
cluster's metallicity is on the lower end of its range then 
this would indicate all elements have been enhanced in this 
cluster. 



JaFu 2 may indicate a depletion of oxygen in NGC 6441. 
The value of log O/H - log O/Hq = -0.93 and the value of 
log Ar/H — log Ar/HQ — —0.72. Since this is a lower metal- 
licity cluster we would expect oxygen as an alpha element 
to be enhanced however it appears to be less enhanced rel- 
ative to argon. This would be consistent with the enhanced 
material having been processed by CNO cycling. 

As GCs age they leak stars into the field via collisions 
and also by tidal stripping. In fact it is estimated that 4- 
9% of halo stars are 2nd generation globular cluster stars 
(jVesperini et al.ll2010l ). If second generation halo stars be- 
come halo stars they could produce PN similar to JaFu 1 
and JaFu 2 and may be GJJC-1. However, this is a small 
fraction of the halo stars and it is quite possible no such star 
has been observed. 

Other places to look for similar PN might be the system 
of satellite galaxies. The GCs u Centauri, M54 and M22 
are all possible captured satellite galaxies and all of them 
have multiple populations (see Piotto 2009). This suggests 
multiple populations might be common in satellite galaxies 
and this model predicts there may be PNe in these satellite 
galaxies similar to the GC PNe. 



5 CONCLUSIONS 

This paper presents a series of models for the expected TP- 
AGB stars in globular clusters. The results of these models 
are compared to the observed abundances of the globular 
cluster PNe and the measured masses of WDs. 

(i) These models suggest the typical unenriched ZAMS 
star at the turn-off will not produce a visible PN, which 
explains small number of PN in GC. The CSPN mass is too 
small to produce a visible PN during the transition from the 
AGB to the white dwarf phase. 

(ii) The PNe JaFu 1 and JaFu 2 are consistent with being 
produced by second generation stars in the globular cluster. 
Both have He/H ratios which are consistent with Fzams ~ 
0.30 which is consistent with results for second generation 
stars from fitting colour-magnitude diagrams to the cluster. 

(iii) Due to its large core mass (0.56-0.58 Mq) it is quite 
possible the PN GJJC-1 is consistent with a third generation 
{Y ~ 0.40) progenitor. 

(iv) The nebula and central star Ps 1 and K648 respec- 
tively may have been produced from a blue staggler. How- 
ever, this interpretation is in doubt since the nebular mass 
of Psl is too small to be explained by this model. The high 
CSPN mass favors this hypothesis but the low nebular mass 
does not favor it. Additional mass-loss due to binary inter- 
action could explain this. 

(v) This model suggests that there should be, in globu- 
lar clusters with different populations with different helium 
abundances, two white dwarf cooling tracks with the sec- 
ond generation star with higher Y producing a slightly more 
massive WD then the first generation star with a lower Y. 
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